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1. Introduction 

Event shapes are among the most powerful observables used for testing QCD [|[, providing, 
for example, some of the best measurements of Og to date Q. However these measurements 
are complicated by the fact that event shapes are unusually sensitive to soft, long-distance 
effects which are expected to be beyond the reach of perturbation theory. These effects give 
contributions to the mean values of the shape variables, so-called power corrections, that 
only fall off like the first power of the centre of mass energy Q, and remain sizable even at 
relatively high energies. Indeed state-of-the-art next-to-leading order (NLO) MS—scheme 
perturbative calculations give a totally unsatisfactory description of the data even for Q 
as high as Mz without adding an additional power-suppressed contribution. 

Unfortunately there is as yet no rigorous framework for computing QCD scattering 
amplitudes non-perturbatively, so a phenomenological model is required to handle the soft 
effects. This can then be combined with NLO perturbative calculations to give a successful 
theoretical description of the data. 

There are different approaches to modelling the power corrections. One possibility 
is to use Monte Carlo event generators to calculate the difference between “parton-level” 
and “hadron-level” distributions of event shapes. More recently an analytic ansatz has 
been proposed, relating the power corrections for many event shape means to a single 
non-perturbative parameter, ao(/^/) |@, @]- This can be interpreted as the mean value of a 
hypothetical universal, infrared-finite coupling aiR(/i) over the range 0 < < fij. Typically 

m is taken to be 2GeV, and the ao values required to fit the data are scattered around 
0.5, with a spread of around 20% (see e.g. the review Ref. [Q| and references therein). 
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Impressively, this holds for event shape means and distributions measured both in e'’'e 
annihilation and deep-inelastic scattering (DIS). 

However, in a recent analysis of e'’'e“ event shape means [^, the DELPHI collaboration 
found that when using an “optimised” choice of renormalization scale the NLO perturbative 
QCD calculation could describe the data very well without any non-perturbative correc¬ 
tions. The idea behind an “optimised” scale is to use some sort of theoretically motivated 
criterion for choosing an observable-dependent renormalization scale ^ (and in general 
scheme) based on features of the perturbative calculations for that observable. The hope 
is that this minimises the effect of higher-order corrections better than the standard pro¬ 
cedure of setting the MS-scheme renormalization scheme to Q. The particular scale choice 
used in the DELPHI analysis was originally proposed by Grunberg; it is the consequence 
of applying the Method of Effective Charges (ECH) |p at NLO. Several different theoret¬ 
ical arguments have been proposed to motivate this scale choice (going under the names 
of the Fastest Apparent Convergence criterion, Renormalization-Scheme Invariant Pertur¬ 
bation Theory [^, Complete Renormalization Group Improvement and Lambda-based 
Perturbation Theory |^]). 

This surprising result motivates investigating the application of “optimised” scales to 
other observables which usually require large power corrections to fit the data. In Ref. ||l^ , 
the Method of Effective charges was applied to the distributions of thrust and heavy- 
jet mass in e'’'e“ annihilation, including a resummation of large logs to next-to-leading 
accuracy. The optimisation did indeed reduce the need for power corrections somewhat, 
although surprisingly the resummation caused problems in the 2-jet region. 


Another class of related observables, which lack the complication of large logarithms, 
are the event shape means defined in the Breit frame of ep DIS. In this case, however, 
the appearance of the proton in the initial state leads to a dependence of the event shape 
means on the proton parton distribution functions (PDFs) and the accompanying factori¬ 
sation scale, M. Ideally we would like to use some optimisation criterion to choose all the 
unphysical parameters, but none of the above-mentioned theoretical arguments simultane¬ 
ously select unique values for both p and M. However, there is an alternative criterion for 
choosing the renormalization scale, the Principle of Minimal Sensitivity (PMS) ||^, that 
does generalise very easily to the factorisation case [1^, IJ, 0 (for recent applications to 
hadron-hadron interactions see Refs. il,IIl)- Moreover, the PMS tends to agree very well 
with the ECH as used in the DELPHI analysis when the latter can be applied. Therefore, 
the purpose of the present paper is to apply the PMS to the choice of renormalization and 
factorisation scales in DIS event shape means, and compare NLO perturbation theory to 
data from the HI collaboration |l^ to see how the required power corrections are affected. 

The plan of this paper is as follows. Section contains definitions of the observables. In 
Section we review the way that dependence on p and M arises in NLO DIS calculations, 
and define the PMS approximations we are going to study. Then, in Section Q we describe 
the methods used in the calculation. Section ^ runs through the analysis in detail for one 
specific observable, the current jet thrust Tc- Results for all the observables discussed in 
Ref. are presented in Section p. Section contains our conclusions. 
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2. Definition of the Observables 


This section briefly summarises the observables which will be considerd in the paper. For 
discussions of cuts and more details about the definition of each observable, see Ref. jl^]. 

In all the following definitions, the momentum of the ith particle is pi and Breit frame 
quantities are indicated by an asterisk. The current hemisphere (CH) is defined as that for 
which 2 * < 0, and the remnant hemisphere (RH) is that for which z* > 0. 

First, we consider five observables defined using particles in the CH only. To ensure 
infrared safety the total energy in the CH (evaluated in the Breit frame) must exceed some 
threshold, taken to be Q/10, for an event to contribute to one of these variables. The 


z*-axis thrust is 


Tr. = l- 


Ei \P*i-n\ 

Ei \Pi*\ 


where n is a unit vector in the z* direction. Similarly, the z*-axis jet broadening is 


Bp = 


The thrust-axis thrust is 


Tp = 1 — max, 


Ei \fi X n| 

E^\m ■ 
E* \n-A 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


E^\m 

where n is now the unit vector with respect to which the maximisation is performed (the 
value which leads to the maximum defines the thrust axis). The C-parameter is defined by 
making use of the tensor 

EiPljPlkl\'Pi*\ 




Ei \Pi 


(2.4) 


with eigenvalues Ai, A 2 and A 3 : 


Cp — 3(AiA2 + A 2 A 3 -|- A 3 A 1 ). 


(2.5) 


Lastly we have the jet mass 


PE 


iErP? 

^E^Er 


( 2 . 6 ) 


Because the jet mass involves the difference between energies and 3-momenta it is especially 
sensitive to hadron mass effects which can lead to additional power corrections beyond the 
scope of the model of Ref. ii- These can be removed by defining a related observable 
in the so-called “E-scheme” |]T] 


(ESch) 

Pe 


E^{El 




^E^Et 


(2.7) 


This is just Pe with each particle replaced by a massless particle of the same energy. Since 
our NLO perturbative calculations are carried out with the approximation of massless 
quarks, the perturbative predictions for will be identical to those for pp'-, only the 

non-perturbative contributions will differ. For the purposes of this paper, values have been 
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obtained for by applying a correction factor from PYTHIA 6.2 to the pE values 

given in Ref. |15]. 

The remaining two variables make use of both the CH and the RH. They are defined 
based on two jet clustering algorithms, the factorisable Jade algorithm and the modified 
kx algorithm. These clustering algorithms allow one to gather the final-state momenta into 
a certain number of current jets, plus the remnant jet, according to a resolution parameter 
y which determines how aggressively momenta are combined. This is done by examining 
the momenta repeatedly, combining those that are nearest according to a distance measure 
specific to the clustering algorithm, until the nearest pair have a separation yij > y. Event 
shapes can be defined from these algorithms by using the value of y at which the event 
goes from being a 2-|-l jet event (two current jets, one remnant jet) to a U-l jet event. 

In the DIS case the distance measures need to be generalised from the e'’"e“ case to 
include a distance between a jet and the remnant. For the factorisable Jade algorithm 


2 F*F;(1 -cos0 ^) 

= Q2 

(2.8) 

2E*XBE;{l-coset) 

y- = g2 

(2.9) 


where i,j indicate jets, r is the remnant and p is the proton. In the same notation the kt 
measures are 


2min(Ef,Ef)(l-cos0*,.) 
= -gl- 

2Ef(l-cos0p 


( 2 . 10 ) 

( 2 . 11 ) 


The related event shapes will be denoted by yjj and yu^ respectively. 


3. Applying The Principle of Minimal Sensitivity to DIS 

All perturbative predictions in QCD are plagued with a dependence on the arbitrary renor¬ 
malization scheme (RS) and scale {p), together referred to as the renormalization prescrip¬ 
tion (RP) 1^. For a process with hadrons in the initial state, predictions depend also 
on the factorisation scheme (FS) and scale (M), which we jointly refer to as the FP. This 
occurs despite the fact that the exact predictions of the theory are completely indepen¬ 
dent of these choices, because truncating the perturbation series prevents cancellations 
between different orders of perturbation theory which are necessary if the dependencies are 
to vanish. 

Because the renormalization and factorisation prescriptions are in principle completely 
arbitrary, they may be varied at will. In the general case this would allow one to obtain 
any desired perturbative prediction for a given observable. Therefore, some method for 
selecting renormalization and factorisation prescriptions is required that will give “reason¬ 
able” quality perturbative predictions. In practice one normally uses the MS RS and FS, 
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and sets ^ = M = Q. In this paper we are going to consider instead picking the scales 


using the PMS |11|. 

The PMS simply states that if we are given an approximation which depends on some 
parameter arising only as part of the approximation procedure (so that the exact result 
is independent of the parameter), the approximation is most believable where it is least 
sensitive to this “unphysical” parameter. This gives a simple way to choose an “optimal” 
value for the parameter (in the ideal case; in general there may be multiple “PMS points” 
and it may be necessary to use some other criterion for choosing between them). 

This principle is not the sort of thing one could imagine proving mathematically; indeed 
it is easy to invent approximations where the PMS fails totally. The sense in which the PMS 
gives an “optimal” approximation is certainly not to be understood in the mathematical 
sense of minimising the difference between the approximate and exact values. Nonetheless 
the PMS can be very successful when applied to real world approximations. 

In perturbative QCD one is in exactly the situation where the PMS is supposed to 
apply. The unphysical parameters label the renormalization and factorisation scales and 
schemes; the full sum of the perturbative series for some observable quantity is independent 
of these parameters but the truncated, partial sum (which is all we have to work with in 
practice) is not. 

The simplest application of the PMS to QCD involves the case of an observable without 
hadrons in the initial state (and without identihed hadrons in the final state). The only 
unphysical parameters are then those that label the renormalization scale and scheme. At 
NLO it suffices to just consider the renormalization scale n, as a change in scheme can be 
absorbed into a rescaling of n- For an observable involving hadrons in the initial state, 
we also have to take into account its unphysical dependence on the parameters that label 
the factorisation scale (M) and scheme. In contrast with the renormalization case, the 
factorisation scheme dependence cannot be absorbed into the scale even at NLO. However, 
it is difficult to formulate a PMS condition for the factorisation scheme in x-space 
so here we will simply neglect the scheme dependence and work in the MS factorisation 
scheme at all times. Hopefully this partial optimisation still represents an improvement 
over the physical scale approach. With this simplification, the relevant parameters are [i 
and M. 

The mean of some DIS observable y depending on the final state X can be expressed 


as 


{v) = 


f y{X)da{ep X,Q) 
f da{ep X, Q) 


(3.1) 


where da{ep X,Q) is the infinitesimal cross-section for the process ep X, and Q is 
the virtuality of the exchanged photon. To compute such a cross-section in perturbation 
theory we need to factorise the process into a soft part, described by the proton PDFs 
fa{^, M), and a hard part, described by a partonic cross-section d{e a ^ X,Q, M) 


da{ep X, Q) 
dX 


/ dUa{i,M) 


da{e a ^ X,Q, M) 
dX ' 


(3.2) 
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As indicated, / and a depend on the unphysical parameter M even before they are ap¬ 
proximated perturbatively; this dependence cancels between them after the integral over 
^ and the sum over a are performed. The dependence of fa{^,M) on M is given by the 
DGLAP evolution equations 


dfa{^,M) 

dlnM 


as{M) 

TT 



PaMz)h{z,M). 


(3.3) 


The Pa^biz) are the splitting functions (for a parton of type b splitting to give a parton 
of type a with momentum fraction z). They are uniquely defined once the factorisation 
scheme is specified. 

To arrive at a NLO approximation, we substitute into Eq. ( p.2| ) PDFs evolving accord¬ 
ing to the NLO splitting functions, and the partonic cross-section expanded to NLO. This 
approximation not only prevents the exact compensation of M dependence between the 
PDFs and the partonic cross-section, but also introduces a dependence on the RP used for 
the expansion. 

In summary, at NLO, our approximation for (y) will have an unphysical dependence 
on both M and y and we can look for a PMS point by requiring 


^(i/)NLO 

dp, 


MPMS 


d{y)KILO 

dM 


= 0 . 


MpMS 


(3.4) 


In general, the stationary point will be a saddle-point in the {p, M) plane. 

This is the most straightforward way to apply the PMS to DIS, but the multiplicity 
of initial states compared to the e'’'e“ case allows for some more involved possibilities. 
In particular, it is interesting to consider the possibility of using different renormalization 
scales in the various partonic channels. This is legitimate as the cross-section for each 
partonic sub-process is separately RP independent, as can be seen by noting that every 
term in the sum over a in Eq. ( |3.2| ) is independent of p (and all other RS parameters). 
However, it must be borne in mind that these cross-sections do depend on the factorization 
prescription - e.g. the dependence on M only disappears after the sum over a is performed. 
Therefore, only this sum is actually physically observable. Now, the usual advice is to 
avoid adding together observables prior to optimisation (this makes sense as the PMS 
point could otherwise arise as a consequence of cancellations between the two observables, 
which is presumably not a good indication of reliability), but these components are not 
observable, so it might not be sensible to apply this rule here. 

The possibility of using different scales for the q'y* and 57 * sub-processes was also 
suggested in Ref. |2^, where the 2 -|- 1-jet cross-section in DIS was studied using a number 
of “optimisation” methods. There the contributions to the jet cross-section for each sub¬ 
process were thought of as being in principle observable (assuming one could distinguish 
quark and gluon jets), because of the characteristically different hnal states: at Born level, 
the quark-initiated “QCD compton” process leads to a qg final state whereas the gluon- 
initiated “boson-gluon fusion” process leads to a qq final state. However, this appears 
to contradict the fact that the cross-sections for these processes depend on the FP. The 
problem arises from considering the final state of the hard process in isolation, when the 
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separation into hard and soft processes is arbitrary and depends on the factorization pre¬ 
scription. For example, as one increases M a part of the contribution from g^* qq where 
the anti-quark is nearly collinear to the proton direction is reinterpreted as a q'j* q pro¬ 
cess with a modified fq{^), obtained by solving the DGLAP evolution equations. Of course 
in the case of a 2 -|- 1-jet cross section with ycut sufficiently large that widely separated 
jets are obtained these ambiguities are lessened, because this implies a large separation 
between the products of the hard and soft parts of the process. 

The reason it becomes interesting here to consider this two-scale PMS procedure is 
that, for M of order the hard scale Q, the quark and gluon channels give quite different 


contributions to the weighted integrals in Eq. (^), to the extent that the NLO contri¬ 
butions are often of different signs. There can be large cancellations between them, and 
it seems plausible that we could get improved results by “optimising” them separately 
(because the cancellations may well be specific to NLO). 

Therefore, we will consider two variants of the PMS in this paper: PMSi where we 
optimise with respect to a single renormalization scale g and the factorization scale M, 
and PMS 2 where we optimise with respect to renormalization scales jiq and gg for the 
quarks and gluons respectively, along with M. It is also possible to assign different scales 
to the different flavours of quarks, but this makes no significant difference to the results 
because the optimisation procedure always chooses very similar scales for them. Therefore, 
for simplicity we confine ourselves to considering g,q and gig. 

For PMSi we can define {y) in the conventional way, by expanding the RHS of Eq. ( |3.1| ) 
in and truncating it at NLO. However, this won’t work for PMS 2 because there are 

two different couplings, and a^{gg), appearing in both the numerator and the 

denominator. Instead we can perform a double expansion in ct^{gq) and a^{fj,g), keeping 
all terms with two or fewer a’s. (We could also truncate the numerator and denominator to 
NLO and use the quotient as our NLO approximation; this makes no qualitative difference 
to the results of this paper). 

In performing the optimisation some very low MS scheme scales can be obtained (even 
down to < IGeV), which raises the question as to whether mass thresholds should be taken 
into account (as well as the more general question of the validity of the results, which we 
hope to assess later by comparison to the data). The problem is that doing so conflicts 
with the idea that /x is a strictly unphysical parameter - indeed, adding thresholds at any 
given quark masses would break the scheme invariance of the PMS results. Another way of 
looking at this is to point out that the optimisation is really being performed with respect 
to the dimensionless parameter r = 61n(/x/A) |]^; varying g at fixed A is just a convenient 
way to achieve this. Therefore, in this paper we treat the number of flavours as fixed to 5, 
as suggested by the value of the physical scale Q. 


4. Calculational Methods 


NLO perturbative predictions for our observables were obtained from the matrix element 


integration Monte-Garlo program DISENT 0.1 |^|, using the interface library NLOLIB 
1^]. Except where noted, the MRST2001E PDF set [|^] was used, and a-^{Mz) was 
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fixed to 0.119 to be consistent with these PDFs. To allow for optimisation with respect 
to M, many runs of DISENT were carried out at intervals AM = IGeV. For every value 
of Q and M, 10® events were generated to ensure the Monte-Carlo integration errors were 
negligible. The /r dependence was implemented by having DISENT compute the coefficients 
of each power of with fi = Q] the /r logs as well as the factors of could then be 
added in later. This means that only one run of DISENT was required for all values of /x. 
However, to allow this to work it was necessary for all the events to be generated at the 
same value of Q, so the mean value of Q was used for each bin rather than integrating over 
the entire width of the bin. The effect of this was estimated by comparing the two results 
for fi = M = Q, giving a correction factor (generally of order 1%) which could be applied 
to the data. 

The power corrections can simply be added to the perturbative prediction for the shape 
means 

{y) = {y)p + (y)pc- (4.i) 

We will present fits both for a simple C\/Q or C 2 lQ^ term, and for the more sophisticated 
universal formula derived in Refs. |^, |^, According to this model 

o.p-i{yi) ~ ^ 

(log(^//x/) K/f3o + l/p)a|jg(/x)] . (4.2) 


(y)pc = 


32 M 


Svr^ 


P 


The overall constant Uy and the integer p depend on the observable. ap-i{pi) is the 
{p — I)th moment of the coupling for 0 < p < pj. A4 ~ 1.49 is the so-called Milan factor 
1^, and K is 


K = 


67 


TT 

T 


-Nf 


(4.3) 


The terms in (y)pc involving are intended to subtract the part of the soft-gluon effects 
that is already included in the perturbative calculations. These are to be evaluated with 
Nf = 5 and with p the renormalization scale used in the perturbative part. Note that 
for our PMS optimised results this will in general differ from Q. In the case of PMS 2 
the appropriate scale is pq as the power corrections do not receive contributions from the 
gluon-initiated sub-process [p^ . 

For all our observables apart from ykf, P = 1 ap has been calculated, so we can 
use this formula to perform a fit for ao- For p = 2 and ap is unknown, which prevents 
a reliable extraction of ai. The coefficients are: 


37r 1 

“tc = 1 o^Cp = o,pE — 2 


(4.4) 


^VfJ 


= 1 Q aBp = 


TT 


-= + R|_o.3069 + 0(1) m 

3“CMW(e-3/4Q) 


(4.5) 


In performing these fits we need to consider both experimental and theoretical sources 
of error. The experimental errors consist of both a statistical and a systematic component; 
lacking knowledge of the proper correlation between the systematic errors we have treated 


- 8 - 









them as uncorrelated, and simply added them to the statistical errors in quadrature, to 
arrive at a composite experimental error. Where the errors are asymmetric, the maximum 
value was taken. 

The MRST2001E PDFs allow a PDF-related error to be estimated by sampling from 
an ensemble of different distributions. Ideally we would repeat the entire analysis with the 
different PDFs, but this is impractical because of needing to rerun DISFNT for each value 
of M. Instead we have estimated the error by comparing the analyses at M = Q only. 
However, the results are so stable with respect to changes of M that this should not be 
too much of a restriction (and in any case, the errors due to uncertainties on the PDFs are 
small). 


5. Case Study: Tc 

In this section we show the results of our analysis, along the lines described in the preceding 
sections, when applied to one typical observable, 1 — Tc = Tc- A summary of the results for 
all the observables is given in the next section. 

In order to arrive at our PMSi predictions, we need to examine the dependence of 
(tc) on ^ and M. This is illustrated in Fig. |^. Note that the PMS points, defined by 
Fq. (|3.4|), are saddle points. Specifically, they are maxima in the /r-direction and minima 
in the M-direction. The actual PMS scales are given in Table 1. 

The effect of choosing these scales over the standard choice // = M = Q is illustrated 
in Fig. 1^. Evidently a substantial power correction is required to fit the data even when 
using the PMS scales. 


Q/GeV 

7.46 

8.8 

14.95 

17.73 

23.75 

36.69 

57.61 

80.76 

/x/GeV 

2.11 

2.45 

4.07 

4.96 

6.98 

12.3 

19.8 

22.4 

M/GeV 

7.87 

7.01 

6.69 

6.25 

5.55 

3.82 

2.72 

2.35 


Table 1: PMSi scales for (tc). 


To evaluate the PMS 2 approximation requires one to look at the dependence of (tc) 
on ^q, fig and M. Some typical examples are illustrated in Fig. ^ The PMS points here 
are maxima in the fiq and fig directions, and may be either minima or maxima in the M 
direction; the relevant scales are listed in Table 2. There is a clear difference between the 
values of fiq and fig, the former are very small and the latter very large. The result of using 
these PMS 2 scales is shown as the dotted line on Fig. |2|. Clearly, at least for large Q, the 
size of the required power correction is substantially reduced, but at low Q the prediction 
now actually overshoots the measurement. It is perhaps not surprising that we run into 
trouble here, as fiq is extremely low, in fact around 0.5GeV~ 2Aj^, so ctjj^{fiq) ~ 0.7. 

To see why PMS 2 makes such a large difference compared to PMSi h is helpful to fix 
M = Q - (it turns out that the effect of optimising M is rather small, and the relation¬ 
ship between PMSi and PMS 2 is simplest when they have a common M). Consider the 
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Figure 1: Dependence of on /x and M for various values of Q. In each case, the PMS 

point is marked. 


coefficients in the perturbative expansion for the integral 


= Aqas{iiq) + + /3olog + 

Agasifig) + [Bg + (5Q\og[-^\Ag \ al{^g). (5.1) 


The coefficients A and B depend on both M and Q. Taking Q = 7.46GeV where the 
difference between PMSi and PMS 2 is most obvious (and fixing M = Q as noted above) 
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Figure 2: Comparison of NLO predictions to HI data for (tc) (combined statistical and 
systematic error bars are shown). The solid curve uses ^ = M = Q. The dashed curve uses the 
PMSi scale choices given in Table 1. The dotted curve uses the PMS 2 scale choices shown in Table 

2 . 


we have 

Ag = 5.29, Ag = 1.44, Bg = 15.4, Bg = -10.3. (5.2) 

On the other hand, if ^g and fig are identified as in PMSi, the coefficients in the pertur¬ 
bation series are simply 


A = Aa + Aa = 6.7A B = Ba + Ba = 5A0. 


(5.3) 


Because of the cancellation in the NLO coefficient B, PMSi sees a series which appears 
to be much more convergent than that seen by PMS 2 , and this lessens the effect of the 
optimisation (this can be understood most easily by recalling the similarity between PMS 
and the Method of Effective Charges, which hxes the scale so that B = 0). Therefore, it 
may be that PMSi underestimates the size of higher orders in the fi = M = Q series. If 
this cancellation does not persist to higher orders, then one would expect PMS 2 to give a 
more realistic estimate of the higher order terms. 

This explanation for the values of the PMS 2 scales is an over-simplification, because 
we aren’t actually optimising the weighted integral Eq. (^), but rather the ratio Eq. (^.Ij). 
However, the total cross-section is convergent enough that these simple considerations do 
capture the essential reason behind the PMS scales shown in Tables 1 and 2. For example. 
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the Effective Charge scales /tec = exp(—i?/(/3o^))Q corresponding to these coefficients 
are 

/i = 4.0GeV (5.4) 

fiq = 0.68GeV fj.g = 2.6 • lO^GeV (5.5) 

which are all rather close to the corresponding PMS scales (although the agreement here 
is perhaps deceptively good and worsens at higher energies, the EC fig actually falling as 
Q is increased). 




Figure 3: Dependence of on ^q, /ig and M for the extreme Q values. The left hand plot 

shows the dependence on fj,q and with M fixed to its PMS value. The right hand plot shows the 
dependence on M and ^q with fig fixed to its PMS value. The PMS points are marked, and the 
numerical values of the PMS scales are listed in Table 2. 


Tc is expected to receive 1/Q power corrections, which we can try to describe either by 
simply adding a term CijQ to the perturbative predictions, or by using Eq. (4.2) which 
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g/GeV 

7.46 

8.8 

14.95 

17.73 

23.75 

36.69 

57.61 

80.76 

/Xg/GeV 

0.50 

0.59 

0.90 

1.1 

1.4 

2.0 

3.2 

4.2 

f^g/GeV 

2.3-103 

6.0-103 

6.5-10'^ 

1.5-10® 

4.1-10® 

1.0-10® 

1.9-10® 

2.2-10® 

M/GeV 

4.95 

5.20 

6.23 

6.41 

6.76 

6.70 

5.79 

4.50 


Table 2: PMS 2 scales for (tc). 


relates the corrections to oq- Because of the rise in the PMS 2 predictions for Q < lOGeV, 
these 1/Q corrections alone cannot compensate for the discrepancy between theory and 
data in this case. In addition, there must be large higher-order ©(a^) effects and/or sub¬ 
leading power corrections oc l/Q^. To compare the size of the power corrections required 
by the different perturbative predictions, it makes sense to exclude low Q data if this gives 
an unacceptably bad fit. This is especially true in view of the fact that the e+e“ data 
examined in Ref. had {Q) > 45GeV, so these additional effects might be important at 
low Q in the e^e“ case also. Therefore, we performed minimum-fits, adding in data 
points from the highest Q downwards until the fit probability fell below 5%. Experimental 
errors were estimated by allowing to vary within 1 of its minimum value. Fitting in this 
way for Ci gives 


^ = M = Q: Gi = 1.23(ll)(20)GeV, Q > 30GeV 
PMSi : Cl = 0.65(2)(2)GeV, Q > 8GeV 

PMS 2 : Cl = 0.18(3)(5)GeV, Q > 14GeV. 

Here and throughout, the first number in brackets indicates the error due to experimental 
and PDF uncertainties; the second number gives an indication of the sensitivity to the fit 
range by showing the size of the shift induced by excluding the lowermost bin. As expected, 
the required power correction is largest for the “physical scale”, and smallest for PMS 2 . 
Of the three predictions, PMSi gives the best description of the data, fitting the largest 
range of data with the best stability. 

Because one expects also sub-leading power corrections to be present, it is interesting 
to introduce e.g. a C2lQ‘^ term into the fit, to see how this affects the conclusions: 

^ = M = Q-. Ci = 1.09(5)(6)GeV, C2 = -4.3(5)(9)GeV2, Q > 7GeV 

PMSi : Cl = 0.82(5)(l)GeV, C 2 = -2.5(5)(2)GeV2, Q > 7GeV 

PMS 2 : Cl = 0.49(5)(0)GeV, C2 = -5.3(5)(l)GeV2, Q > 7GeV 

where the errors are strongly correlated. Unsurprisingly, the C2lQ‘^ term allows even the 

low Q data to be correctly described. The basic fact that the optimisation reduces the 
need for 1/Q power corrections does seem to survive. 

Fitting for oq gives 

^ = M = Q: ao = 0.524(7)(5), Q > 8GeV 
PMSi : ao = 0.596(6)(7), Q > 7GeV 
PMS 2 : ao = 0.614(8)(16), Q > 14GeV. 
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It may perhaps seem surprising that the oq values are larger for the PMS scale choices, 
even though the size of the power corrections seems to be reduced, as is reflected in the 
Cl values. This is a consequence of using the PMS scales in Eq. (4.2), which increases the 
perturbative contribution and requires larger oq to compensate. 

Finally, we can examine the effect of assuming a different value for For consis¬ 

tency, this requires the PDFs to be changed. Fig. Q shows the effects on (tc) of choosing 
ay^{Mz) = 0.1125 and 0.1225, using PDFs from Ref. [|^]. To save computer time M was 
fixed to Q throughout. This makes little difference to the results - indeed, this can now 
easily be seen by comparing the central graph of Fig. ^ with Fig. which only differ in 
that M is fixed to Q in the former and optimised in the latter. 

The results of performing power corrections fits for these values of are summarised 
in Table 3. The effect of the optimisation is still to reduce the size of the power corrections 
(although less substantially for low as would be expected). 


Ci/GeV 


/i Q 

PMSi 

PMS 2 

0.1125 

1.43(11)(27) 

> 30 

0.74(2)(2) 

> 8 

0.51(4)(6) 

> 16 

0.119 

1.23(11)(20) 

> 30 

0.53(2)(3) 

> 7 

0.13(3)(5) 

> 14 

0.1225 

0.87(6)(24) 

> 20 

0.44(2)(2) 

> 7 

-0.09(2)(1) 

> 8 

ao 

0.1125 

0.661(33)(69) 

> 30 

0.586(7)(5) 

> 8 

0.559(6)(4) 

> 7 

0.119 

0.524(7)(5) 

> 8 

0.596(6)(6) 

> 7 

0.575(8)(15) 

> 14 

0.1225 

0.511(6)(9) 

> 7 

0.615(6)(2) 

> 7 

0.519(18)(5) 

> 20 


Table 3: Fits to (tc) for Ci and cxq, with various value of aj^{Mz). Note that M = Q throughout. 
The first number gives the best fit value and numbers in brackets indicate errors in the last digits 
(due to experimental and PDF uncertainties) and the variation arising from removing one bin from 
the fit. The notation > Q indicates the range of Q that could be fitted before the indicated a 
fit probability of < 5%. 


6. Results 

In this section we summarise results for all the observables studied in Ref. |^. The 
perturbative predictions are compared to data in Figs. ^ and Tc and Cp show similar 
features: the PMSi predictions are somewhat closer to the data than the physical scale ones, 
and the PMS 2 predictions are a lot closer, except they become too large at low energies. 
The PMS predictions for rp are closer to the data than those with the /r = M = Q, with no 
evident breakdown at low Q. Bp and pp also have reasonable low Q behaviour but with a 
lesser improvement of the fit. Turning to though, the improvement (especially for 

PMS 2 ) is much more substantial. The jet transition parameters, yfj and move away 
from the data when the scales are optimised. 

Table 4 shows fits for a Ci/Q power correction for all observables except and a 
fit for a power correction for y^^ (y^^ is the only observable whose leading power 
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Figure 4: Effect of varying on the values of (tc). The top graph has aj^[Mz) = 0.1125, the 
middle has ajj^{Mz) = 0.119, and the bottom has a-^{Mz) = 0.1225. Throughout we have fixed 
M = Q. The solid line has fj, = Q, the dashed line has /i picked by PMSi, and the dotted line has 
/j, picked by PMS 2 . Comparing the central graph with Fig. ^ shows the effect of optimising M. 


correction is expected to be oc l/Q^). Table 5 shows the results of fitting for a power 
correction based on Eq. (^) for all observables bar These fits allow us to see to what 
extent the discrepancy between the perturbative predictions and the data visible in Figs. ^ 
and P can actually be described as a power correction. As was noted in Section PMS 2 
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Figure 5: Predictions of our three perturbative predictions compared to data. The solid line uses 
the scale choice ^ = M = Q, the dashed curve uses PMSi and the dotted curve PMS 2 . Data are 
shown from Ref. |Q . 


fails to give a good description of Tc at low Q. This is also the case for Cp and pp (and 
^iESch)y Although not obvious from Figs. ^ and the “physical scale” predictions don’t 
describe the low energy data for these observables very well either. In fact, the best overall 
fits seem to be those that use PMSi. The other observable with a large discrepancy at 
low Q, T/fcj, can be described quite well by any of the scale choices provided we add the 
expected 1 /Q^ power correction. 


Ci/GeV 

Obs 

p = M = Q 

PMSi 

PMSa 

Tp 

0.73(7)(5) 

> 7 

0.45(7)(2) 

> 7 

0.33(7)(1) 

> 7 

Tc 

1.23(11)(20) 

> 30 

0.65(2)(2) 

> 8 

0.18(3)(5) 

> 14 

Bp 

0.79(3)(2) 

> 7 

0.50(3)(7) 

> 7 

0.37(3)(8) 

> 7 

Cp 

4.17(29)(92) 

> 30 

2.26(7)(15) 

> 14 

0.76(49)(24) 

> 16 

PE 

0.92(6)(13) 

> 30 

0.67(4)(12) 

> 20 

0.49(4)(13) 

> 20 

VfJ 

-0.11(3)(3) 

> 7 

-0.11(3)(2) 

> 7 

-0.18(3)(3) 

> 7 

(ESch) 

Pe 

0.65(5)(12) 

> 30 

0.38(2)(2) 

> 14 

0.21(2)(4) 

> 16 

Ca/GeV^ 

Vkt 

-2.70(60)(87) 

> 7 

-6.03(60)(24) 

> 7 

-6.30(60)(6) 

> 7 


Table 4: Fits for Ci and C 2 to all observables. See Table 3 for an explanation. 
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Figure 6: Predictions of our three perturbative predictions compared to data. The solid line uses 
the scale choice ^ = M = Q, the dashed curve uses PMSi and the dotted curve PMS 2 . Data are 
shown from Ref. |Q . 


ao 

Obs 

^ = m = q 

PMSi 

PMS 2 

Tp 

0.519(22)(2) 

> 7 

0.533(22)(20) 

> 7 

0.559(22)(27) 

> 7 

Tc 

0.524(7)(5) 

> 8 

0.596(6)(7) 

> 7 

0.614(8)(16) 

> 14 

Bp 

0.568(17)(20) 

> 8 

0.447(13)(9) 

> 7 

0.443(13)(6) 

> 7 

Cp 

0.465(7)(12) 

> 16 

0.527(3)(6) 

> 7 

0.557(7)(23) 

> 16 

PE 

0.808(38)(66) 

> 30 

0.715(10)(15) 

> 14 

0.736(7)(9) 

> 7 

VfJ 

0.264(18)(22) 

> 7 

0.261(10)(17) 

> 7 

0.278(12)(19) 

> 8 

—(ESThT 

Pe 

0.642(34)(56) 

> 30 

0.572(7)(15) 

> 8 

0.611(7)(10) 

> 8 


Table 5: Fits for ao- See Table 3 for an explanation. 


7. Conclusions 


In this paper we have studied how power correction fits to event shape means in DIS are 
affected by choosing the factorization and renormalization scales according to the Principle 
of Minimal Sensitivity. In doing this, two different prescriptions were adopted: PMSi, 
where the unphysical parameters were taken to be and M and PMS 2 , where different 
values of were used for the quark- and gluon-initiated sub-processes. The motivation 
behind PMS 2 was to avoid underestimating the effect of higher order corrections because 


of the cancellations between the q'y* and 57 * sub-processes at NLO (illustrated in Eq. (5.2)). 

PMSi gives results that are pretty close to those found using the conventional choice 
fj, = M = Q. However, it does improve the quality of the power correction fits (see Tables 
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4 and 5). PMS 2 gives perturbative results that are substantially closer to the data at high 
values of Q, but which deviate from it at low Q (see Figs. and If we exclude this low 
Q region from the fits (as in Tables 4 and 5), PMS 2 requires smaller power corrections to 
fit the data than does either PMSi or the choice /i = M = Q. 

In carrying out these calculations we have used the MRST2001E PDF set. It is possible 
that problems could arise from using optimisation with PDFs that were obtained without 
optimisation. If, say, switching to the PMS scale has a consistent effect on a number of 
observables, refitting for the PDFs using the PMS would cause them to change so as to 
counteract the effect of the optimisation. Unfortunately, optimising enough observables to 
be able to fit for the PDFs would be a very large undertaking, and is beyond the scope of 
this paper. One might hope that this is not actually an issue, because event shapes have 
particularly large NLO corrections and hence are probably more sensitive to optimisation 
than most other observables. 

The motivation for this study was to determine whether “optimised” scales could 
significantly reduce the need for power corrections to DIS event shape means as they 
appeared to do for their e'''e“ counterparts [^. This would suggest that what appears at 
NLO to be a power correction is just the combined effects of higher order perturbative terms 
(although why genuine non-perturbative corrections would be so small is quite mysterious). 
PMSi does not support this idea; PMS 2 does to some extent, but is poorly behaved at low 
energies - however, these energies are much lower than those studied in [^. So if PMS 2 
provides a better estimate of higher order terms in the perturbation series than PMSi, 
it could be that the conclusions of Ref. do extend to DIS event shape means. In this 
case, NNLO effects and/or power corrections would become important in both processes 
at Q ~ 20GeV. 

To try to determine which optimisation is most accurate one can turn to the data. 
Unfortunately, fitting for a power correction largely absorbs the differences between them, 
leaving only the fact that PMSi gives a better description of the low energy data due to 
it lacking the rapid growth of PMS 2 (visible for example on Fig.^). Another way to test 
the different optimisations would be to perform a PMS analysis of the ZEUS data, some of 
which is binned in both Q and xb However, as long as we only have NLO calculations 
to work with it will be difficult to be sure which, if any, optimisations are best (although 
the overall consistency of the Method of Effective Charges analysis in Ref. Q is highly 
suggestive). Once NNLO computations become available for these event shape means it 
will be possible to measure the convergence of the optimised approximants and compare 
them to each other and to the convergence in the MS-scheme. This should provide better 
guidance in choosing the scheme for these observables, and allow us to more thoroughly 
test if the apparent power corrections really can be mimicked by an optimisation of the 
scheme. 
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